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ABSTRACT 


Child  and  Langmuir  introduced  a  solution  to  space 
charge  limited  emission  in  an  infinite  area  planar  diode. 
The  solution  follows  from  starting  with  Poisson''  s  equation, 
and  requires  solving  a  non-linear  differential  equation. 
This  approach  can  also  be  applied  to  cylindrical  and 
spherical  geometries,  but  only  for  one-dimensional  cases. 
By  approaching  the  problem  from  Coulomb' s  law  and  applying 
the  effect  of  an  assumed  charge  distribution,  it  is 
possible  to  solve  for  space  charge  limited  emission  without 
solving  a  non-linear  differential  equation,  and  to  limit 
the  emission  area  to  two-dimensional  geometries.  Using  a 
Mathcad  worksheet  to  evaluate  Coulomb' s  law,  it  is  possible 
to  show  correlation  between  the  solution  derived  by  Child 
and  Langmuir  and  Coulomb's  law. 
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I .  INTRODUCTION 

A.  APPLICATIONS  OF  VACUUM  DIODES 

Recent  emphasis  on  missile  defense  and  the  expansion 
of  electronic  warfare  has  increased  a  need  for  high-energy 
electromagnetic  weapons.  Directed  energy  microwave  weapons 
require  high-energy  electron  beams  to  achieve  the  power 

levels  adequate  for  standoff  utilization,  which  can  reach 

the  gigawatt  range  for  peak  power,  which  is  far  beyond  the 
power  production  of  a  semiconductor  device. 

One  method  of  creating  electron  beams  is  through 
vacuum  diodes.  In  a  vacuum  diode,  an  electron  cloud  is 

created  at  a  cathode  and  accelerated  across  a  vacuum  gap  to 
a  collector  anode.  This  can  be  done  through  thermionic 
emission,  field  effect  emitters,  or  explosive  plasma 
emitters.  For  the  pulse  power  requirements  of  high  power 
microwaves,  explosive  emitters  ensure  that  enough  electrons 
will  be  available  to  create  currents  necessary  for  high 

power  devices  within  the  desired  turn-on  time. 

Experimental  work  on  diode  emitters  is  an  ongoing 
subject,  and  is  focused  on  diode  geometries  and  materials. 
Suppression  of  plasma  effects,  cathode  longevity,  and  the 
use  of  micro-fiber  cathodes  are  all  areas  of  interest  in 
the  creation  of  fieldable  high-energy  weapons.  However, 
the  electromagnetic  effects  of  large  current  and  charge 
densities  in  such  a  diode  are  not  entirely  understood.  One 
effect  is  the  screening  process  from  the  vacuum  charges  of 
the  electric  field  produced  by  the  diode. 

B.  DEFINITION  OF  SPACE  CHARGE -LIMITED  EMISSION 

At  high  current  densities,  the  electric  field  due  to 
the  electrons  and/or  plasma  in  the  diode  cannot  be  ignored. 
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Since  the  operation  of  these  diodes  is  in  the  high  current 
regime,  this  effect  on  the  total  electric  field  is 
important . 

In  a  parallel-plate  vacuum  diode,  the  electric  field 
produced  by  the  diode  is  a  function  of  gap  distance  and  the 
potential  difference  applied.  The  vacuum  electric  field  is 
constant  and  uniform.  In  this  case,  the  current  would  be 
solely  limited  by  the  capacity  of  the  cathode  to  emit 

charge.  However,  experiment  shows  that  the  current 

produced  by  a  real  diode  is  indeed  finite,  regardless  of 
how  much  charge  the  cathode  can  emit.  The  reason  for  this 
is  the  decrease  of  the  electric  field  due  to  the  charge  in 
the  diode  gap. 

As  electrons  enter  the  diode,  the  charge  creates  a 
field  opposite  to  the  field  applied  by  the  diode.  If 

emission  from  the  cathode  is  allowed  to  continue,  the 

electric  field  of  the  charges  can  eventually  cancel  the 

diode  field.  At  this  point,  the  charges  emitted  by  the 
cathode  are  no  longer  accelerated  away  from  the  cathode  and 
current  is  limited.  This  is  space  charge-limited  emission. 
The  current  produced  by  the  diode  is  not  limited  by  the 
ability  of  the  cathode  to  produce  charge,  but  by  the 
canceling  effects  the  current  has  on  the  total  electric 
field  in  the  diode. 

This  limit  has  been  explored  since  1911.  However,  no 
closed  analytic  solution  exists  for  finite  diodes,  which  of 
course,  all  practical  diodes  are.  Since  ultimately,  the 
power  of  an  electromagnetic  weapon  is  limited  by  the  power 
produced  by  the  electron  beam,  this  limit  is  an  important 
factor . 
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II.  CHILD-LANGMUIR  LAW  FOR  SPACE  CHARGE -LIMITED 

EMISSION 


A.  GEOMETRY  OF  AN  INFINITE  VACUUM  DIODE 

The  derivation  of  space  charge-limited  emission  is  a 
corner  stone  of  plasma  physics.  Space  charge-limited 
emission  is  defined  as  driving  the  electric  field  at  a 
cathode  emitter  to  zero,  preventing  any  increase  in  the 
current  density.  The  first  form  of  the  well-known  Child- 
Langmuir  Law  was  published  in  1911.  Beginning  with 
Poisson' s  equation  and  the  one-dimensional  boundary 
conditions  of  a  parallel  infinite  plate  diode,  the  maximum 
current  density  can  be  found. 

The  geometry  of  the  parallel  infinite  plate  diode 
consists  of  an  infinite  grounded  cathode  and  an  anode  held 
at  some  potential,  VQ.  The  two  plates  of  the  diode  are  some 
distance,  D,  apart.  Current  is  then  allowed  to  flow  across 
the  vacuum  gap.  When  the  space  charge-limited  condition  is 
met  at  the  emission  surface,  the  electric  field  produced  by 
the  current  in  the  gap  will  balance  the  electric  field 
produced  by  the  potential  difference  of  the  two  plates. 
This  geometry  is  reflected  in  the  following  diagram. 
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charges  are  being  created  or  destroyed  in  this  steady-state 
problem,  the  relation  is  simply: 

J(z)  =  p(z)*v(z)  (2) 

Then  the  charge  density  is  expressed  as: 


p(z)  = 


v(z) 


(3) 


The  conservation  of  energy  is  used  to  express  the 
velocity  in  terms  of  known  parameters  of  the  diode: 

PE(0)  +  KE(0)  =  PE(z)  +  KE(z )  (4) 


Since  the  charge  is  assumed  to  be  emitted  with  zero 
velocity,  the  kinetic  energy  at  the  cathode  is  zero.  Since 
the  cathode  is  grounded,  the  potential  energy  at  the 
cathode  is  also  zero,  neglecting  gravitation  effects. 
Therefore,  when  expressions  for  the  potential  and  kinetic 
energy  of  an  electron  are  substituted: 


0  =  —e*  O(z)  + 


1 

2 


me  *  v(z)2 


(5) 


Solving  for  velocity  yields: 


v(z)  = 


(6) 


Substitution  of  the  expressions  for  charge  density  and 
velocity  into  the  Poisson  equation  results  in  an  expression 
relating  potential  and  current  density.  The  expression  is 
a  second-order  non-linear  differential  equation. 

d2®(z)  _  J(z)  p7*  1 
8 z2  s0  V  2e  ®(z) 

The  boundary  conditions  that  lead  to  the  space  charge- 
limited  solution  are  the  potentials  at  the  anode  and  the 
cathode,  and  the  condition  that  the  electric  field  at  the 


cathode  be  zero. 

0(0)  =  0 
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(8a) 


®{D)  =  V0 


(8b) 


ao(p) 

dz 


(8c) 


The  last  boundary  condition  is  imposed  because  we  are 
solving  for  the  space  charge-limited  emission  case.  Using 
the  preceding  boundary  conditions,  an  outline  of  the 
solution  to  the  differential  equation  follows.  Multiplying 
both  sides  by  the  derivative  of  the  potential  and 
integrating : 


^  dz  dz 2  J  a0  V  2e  ^  dz 


(9) 


The  integrals  are  of  the  form: 


(  i 


j  (u  *  du )  =  j  u  2  *  du 


(10) 


Solving  these  integrals  yields: 


ao 

dz 


=  -  —  J—*2*®2  +C 
£n  V  2e 


(11) 


Here  the  constant  of  integration  is  zero  from  the 
boundary  condition  that  the  electric  field  is  zero  at  the 
cathode  surface.  Taking  the  square  root  of  both  sides: 


—  =  -2  L 

dz  \  sn 


(  m,  ^ 


1/4 


v2  e  j 


(12) 


Rearranging : 


®  1/450>  =  2* 


f  V/4 
'  m ,  ' 


\  2e  j 


*  dz 


(13) 
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Integrating  this  differential  equation: 


lo3/4 

3 


=  2*  I—  * 


f  \1/4 
'  m ,  ' 


v2  ej 


■■z  +  C 


(14) 


Again,  the  boundary  condition  of  the  grounded  cathode 
forces  the  constant  to  zero.  By  applying  the  anode 
boundary  condition,  one  can  solve  for  J,  the  Child-Langmuir 
Space  Charge  Limit  for  emission  density: 

3/2 


2e  .  K 
111 


D- 


(15) 


This  is  the  maximum  current  density  that  can  be 
produce  from  a  parallel  infinite  plate  diode,  given  a 
potential  difference,  VQ,  and  spacing,  D.  Given  this 
current  density,  substitution  yields  solutions  to  the 
following  quantities: 

Potential : 


0(z)  =  V0 


f  z^n 


\Dj 

Electric  Field: 


(16) 


\E\ 


4  V 

__  * _ o_  ^ 

3  D 


izV/3 


\Dj 


Charge  Density: 


K_ 

D2 


VDJ 


(17) 


(18) 


C.  LIMITATIONS  OF  CHILD-LANGMUIR  DERIVATION 

Although  the  Child-Langmuir  approach  solves  for  the 
space  charge  limit,  it  is  only  valid  for  the  one¬ 
dimensional  case.  It  has  been  generalized  to  infinite 
cylindrical  and  spherical  geometries,  but  it  has  not  been 
solved  for  a  finite  emission  area.  Instead  of  approaching 
the  problem  from  the  Poisson  equation,  which  generates  a 
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non-linear  differential  equation,  if  the  integral 


Coulomb' s 
developed 


law  is  used,  the  space  charge 
by  integration. 


limit 


form  of 
can  be 
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III.  COULOMB ' S  LAW 


A.  COMPARING  POISSON'S  EQUATION  WITH  COULOMB'S  LAW 

In  the  diode  configuration  for  the  space  charge 
limited  formulation,  it  is  assumed  that  the  diode  is  in  the 
steady  state  condition,  and  does  not  vary  with  time.  With 
the  addition  of  a  background  confining  magnetic  field, 
these  conditions  reduce  Maxwell's  Equations  for 
electrodynamics  down  to  Gauss's  Law: 


V-E(r)  =  ^_  (19) 

One  result  of  Gauss's  law  is  the  Poisson  equation, 
from  which  the  Child-Langmuir  solution  for  space  charge- 
limited  flow  is  derived.  As  a  definition,  the  electric 
field  is  the  negative  of  the  gradient  of  the  scalar 
potential : 


E(r)  =  —V®(r)  (20) 

Using  the  definition  of  the  Laplacian  operator: 

V  (-VO)  =  -V2®  (21) 


Substitution  into  Gauss's  Law  gives  the  Poisson 
Equation : 

V2®  (f)  =  -PQ  (22) 

Gauss's  Law  also  generates  the  integral  form  of 
Coulomb's  law.  Substituting  the  integral  of  a  Dirac  delta 
function  into  the  expression  for  p: 

\ — f  4ftS3(r  -  r')p(r')dr'  (23) 

One  form  of  the  delta  function  is: 
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4nS 3  (r  —  r')  =  V 


(r  -  r') 


Substituting  for  the  Dirac  delta  function: 


— - —  [  4 7iSi  (r  -  r')p{r')dr' 
4  nsa  J 

=  1  J 

4ns 0  J 

<1 

Z' 

t  Gt, 
1  1 

-  r’) 

3  p{r')dr 

r')  J 

’  (25) 

Since 

the  integral  in 

the 

above 

equation 

is  with 

respect  to 

r' ,  and  the  divergence  is  with 

respect 

to  r ,  the 

divergence 

can  be  placed  outside 

the 

integral , 

and  the 

constant  placed  inside.  From  Gauss's  law,  this  is  equal  to 
the  divergence  of  the  electric  field: 


V  •  E(r)  =  vj- 


1  (r  -  r') 
K£o  (r-r') 


p(r')dz' 


This  must  be  hold  for  any  r,  so  the  arguments  must  be 
equal.  Thus,  Coulomb's  law  is  derived  from  Gauss's  law: 


E(h  =  \- 


1  ( r  -  r') 

n£o  Cr-r') 


p(r')dz' 


Therefore,  Coulomb's  law  and  Poisson's  equation  are 
two  different,  but  equally  valid  forms  of  Gauss's  law.  The 
results  obtained  by  applying  either  equation  to  the 
conditions  of  the  space  charge-limited  diode  should  be 
equivalent . 

B.  COULOMB'S  LAW  GEOMETRY 

Coulomb' s  Law  in  integral  form  is  an  extension  of 
Coulomb's  force  law  for  two  point  charges.  Evaluating 
Coulomb' s  law  results  in  an  expression  for  the  electric 
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field  produced  by  a  distribution  of  charges.  The  integral 
form  is: 

(28) 

In  this  case,  the  electric  field  can  vary  depending  on 
the  position  where  it  is  measured,  the  field  point,  rf  .  It 
is  also  dependent  on  the  magnitude  and  position  of  the 
charge  distribution.  The  vector  A  is  the  displacement 
between  the  source  point,  rs  ,  and  the  field  point. 

l  =  (r/-0  (29) 

The  geometry  is  reflected  in  the  following  diagram. 


Field 


Figure  2 .  Coulomb  Diode  Geometry 

The  differential  charge  element,  dq,  can  be  expressed 
as  the  product  of  the  charge  density,  p,  and  a  differential 
volume  element.  In  cylindrical  coordinates,  this  becomes: 

dq  =  p(rs )  *  rs  drs  d  0s  dzs  (30) 

To  simplify  the  analysis  of  Coulomb's  law,  the  beam 
can  be  assumed  to  be  confined  horizontally  by  an  external 
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magnetic  field,  thus  require  only  the  analysis  of  the 
vertical  component  of  the  electric  field.  In  this  case: 

X  Z=(zf-zs)z  (31) 

Although  the  directionality  of  X  and  E  have  radial 
and  azimuthal  dependence  removed  by  the  confinement 
assumption,  it  does  not  change  the  geometry  of  the  problem. 

The  magnitude  of  X  remains: 

X  =^r/  +rs2  -2r/rs[sin(^/)sin(^)  +  cos(^/)cos(^)]+(z/  -zj2  (32) 

Now  that  we  have  found  expressions  for  all  the 
components,  the  ability  of  software  such  as  Mathcad  to 
numerically  solve  integrals  becomes  attractive. 
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IV.  MATHCAD  WORKSHEET 


A.  CHILD-LANGMUIR 

Mathcad  is  a  commercially  available  numerical 
calculation  program.  Using  Mathcad' s  symbolic  interface, 
it  is  possible  set  up  a  worksheet  that  compares  the 
electric  field  predicted  by  Child-Langmuir  and  from  the 
integral  form  of  Coulomb's  law.  If  the  two  correspond,  the 
flexibility  of  the  worksheet  will  allow  for  the 
substitution  for  various  charge  distribution  geometries. 

The  first  inputs  into  the  worksheet  are  the  global 
variables  that  control  the  initial  geometry  of  the  problem. 
For  the  initial  purpose  of  comparing  Child-Langmuir  and 
Coulomb's  Law,  the  gap  distance,  D,  and  potential 
difference,  VQ,  are  constant  factors,  and  can  be  set  to 
unity  to  simplify  the  analysis.  The  same  is  true  of  the 
permittivity  of  free  space,  s0. 

The  geometry  of  the  Child-Langmuir  diode  is  two 
infinite  plates.  The  corresponding  parameter  is  an 
infinite  radius,  R. 

For  comparison  the  results  of  the  Child-Langmuir 
solution  are  defined.  The  functional  form  of  the 
potential,  ®(z),  is  given.  Then  the  Child-Langmuir 
electric  field,  Ecl(z),  is  the  one-dimensional  divergence, 
or  the  derivative  with  respect  to  z.  The  vector  component 
of  the  electric  field  in  the  z-direction  is  denoted  by  a 
subscript  z.  Eclz(z)  is  plotted  for  comparison  with  the 
electric  field  calculated  by  Coulomb's  law. 
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Figure  3 . 


Child-Langmuir  Electric  Field 


B.  COULOMB'S  LAW 

The  charge  density  follows  as  the  divergence  of 
Ecl(z)  .  This  charge  density  will  be  used  in  the  Coulomb's 
law  integral.  However,  a  distinction  must  be  made  be  in 
this  transition  between  the  source  point,  zs,  and  the  field 
point,  zf.  In  Child-Langmuir,  the  electric  field  is 
evaluated  at  the  same  z  as  the  charge  density.  For 
Coulomb' s  law,  the  charge  density  is  integrated  over  all 
the  sources  to  give  the  field  at  a  distinct  point. 

The  Coulomb  integral  is  integrated  over  all  0S, 
therefore  the  choice  of  9f  is  arbitrary.  Choosing  0f  to  be 
zero  simplifies  the  expression  for  the  displacement  vector. 
The  magnitude  of  the  displacement  vector  becomes: 


(33) 


As  a  first  check  to  the  model,  the  conditions  of  the 


Child-Langmuir  diode  are  applied.  The  symmetry  of  the 
infinite  diode  allows  rf  to  go  to  zero.  This  further 


simplifies  the  displacement  vector,  removing  azimuthal 
dependence . 


(34) 


The  evaluation  of  the  Coulomb  integral  can  be  plotted 
as  a  function  of  zs. 
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Figure  4 .  Coulomb  Electric  Field 


C.  COMPARING  COULOMB'S  LAW  AND  CHILD -LANGMUIR 


0  z,z  z  ^  D 


Figure  5.  Comparison  of  Child-Langmuir  and  Coulomb' 


Law 
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It  is  evident  from  the  plot  that  the  limits  of  Coulomb 
field  differ  from  the  limits  of  the  Child-Langmuir  field. 

2  V 

Both  the  upper  and  lower  limit  are  too  high  by . 

3  D 

2  V 

Subtracting - from  the  entire  Coulomb  field,  it  is  seen 

3D 

that  the  Coulomb  field  differs  from  the  Child-Langmuir 
field  by  this  constant. 

The  above  plot  shows  that  the  Coulomb' s  law  integral 
and  the  Child-Langmuir  solution  differ  only  by  a  constant. 

D.  ADJUSTING  THE  COULOMB  FIELD 

That  the  Coulomb  integral  differs  should  not  be  a 
surprising  result.  The  integral  only  accounts  for  the 
effect  of  the  charge  in  the  diode  gap.  However,  the  total 
electric  field,  as  predicted  by  Child-Langmuir,  is  also 
influenced  by  the  potential  difference  of  the  conducting 
plates . 

One  approach  to  finding  the  effect  of  the  conducting 
plates  is  the  method  of  images.  However,  the  solving  for 
the  images  of  a  charge  distribution  on  two  conductors  that 
are  held  at  different  potentials  proved  beyond  the  scope  of 
this  thesis. 

Another  approach  is  to  apply  the  boundary  condition  of 
the  potential  difference  across  the  gap. 

K  =  corrected  (Z  f  )dZ  f  (35) 

Here,  Ezcorrected  ( Zf )  is  a  sum  of  the  electric  field  from 
the  space  charge,  and  a  correction  factor,  the  electric 
field  produced  by  the  conducting  cathode  and  anode. 

Ezcorrected  (Z/  )  =  ^ z  (Z /  )  +  ^ zplates  (Z/  )  (36) 

Integrating  the  space  charge  electric  field: 
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(37) 


f  EAzf)dzf=~V0 

Subtracting  from  the  total  potential  difference: 

{0°  E:plates  (zf  )dzf  =  - 1 V0  (38) 

From  this,  the  correction  to  electric  field  from  the 
plates  is: 

E,pa,l.^)  =  ~~  (39) 


This  is  exactly  the  magnitude  of  the  difference  of  the 
Coulomb  field  from  Child-Langmuir .  This  formulation  shows 
a  correspondence  between  the  predictions  of  Child- 
Langmuir'  s  solving  of  a  differential  equation  and 
integration  using  Coulomb's  Law. 

E.  GENERALIZING  THE  COULOMB  INTEGRAL 

This  integral  formulation  has  the  same  limitations  as 
the  differential  equation  solution,  in  that  it  only  holds 
for  an  infinite  area  diode.  To  be  effective  for  solving 
for  a  finite  diode,  the  formulation  must  account  for  off- 
axis  values  of  rf,  introducing  an  azimuthal  dependence. 

The  generalized  Coulomb  integral  does  not  have  simple 
symmetries  to  allow  for  contraction  of  a  dimension  and  is: 


P(Z,Y 


IS  A  H 

•  '  ill 


_ (zf  - z, ) _ 

ylr/+rs2-  2rfP  [cos^ )]  +  (zf  ~  zs ) 
+  rs2  -  2  rfrs  [cos(0, )]  +  (z,  -  z,  )2 1 


rdddrdz  (40) 


In  this  form,  Mathcad  is  unable  to  reliably  solve  for 
Ez(Zf)  over  the  range  0  <  Zf  <  D. 
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V.  CONCLUSIONS  AND  FUTURE  WORK 


The  goal  of  this  thesis  was  to  show  a  correspondence 
between  the  space  charge  limited  conditions  predicted  from 
the  Child-Langmuir  law,  and  the  conditions  predicted  from 
Coulomb's  law.  In  this  respect,  the  Mathcad  worksheet  was 
successful.  However,  it  should  be  noted,  that  assumptions 
imposed  on  the  Coulomb  integral  limit  the  worksheet  to  the 
one-dimensional,  infinite  diode  case. 

The  assumptions  that  limit  the  Mathcad  worksheet  are 
assuming  azimuthal  symmetry  of  the  charge  distribution 
about  the  field  point.  For  a  finite  distribution,  this 
constrains  the  field  point  to  be  on  the  centerline  of  the 
distribution.  This  limitation  does  not  allow  for  the 
exploration  of  the  edge  effects  of  a  finite  beam,  where  the 
maximum  current  density  exceeds  the  Child-Langmuir  limit. 
These  "wings"  could  vastly  affect  the  performance  of 
practical  vacuum  diode,  as  much  more  power  would  be  produce 
than  that  predicted  by  the  simple  infinite  approximation. 

Opportunity  for  future  work  exists  in  extending  the 
worksheet  to  allow  Mathcad  to  solve  for  arbitrary 
geometries.  Including  the  azimuthal  dependence  would 
generalize  the  Coulomb  integral,  making  the  worksheet  a 
powerful  tool  for  exploring  the  electric  field  of  a  vacuum 
diode . 
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APPENDIX 


The  following  pages  are  a  rich  text  copy  of 
Mathcad  worksheet  used. 


the 
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Global 


00  Subscripts  onA,  and  E  represent  vector 

D  =  l 
vo"! 

Eosl 
5=  .001 

Child-Langmuir: 


Potential: 

4 

®(z)  :=Vo(| 

Electric  Field: 


Ecl^z)  :=— — <J>(z) 
d  z 


Ecl^z)  -> 


1 

3 


Charge  Density: 


Coulomb's  Law: 

Infinite  Diode  Case: 


Displacement  Vector: 


%’Vrs)  :=J(rs)2  +  (zf“zs)2 


Charge  Density: 


p(zs)  := 


-4 


Electric  Field: 


Eh)  ■= 


rD 

rR 

fj_\ 

,2&o) 

0  J 

0 

^(zf>zs>rs)2  ) 


drsdzs 


0  dependence  removed  due  to 
azimuthal  symmetry. 


Ez(.0003)  =  0.581 
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^corrected  (?f) '  ^z{7f) 


2- V 

3- D 


Space  Charge  Effect  on  Diode: 
fD 

V:=  e7(t^  zf 

J.oioi 


V=  -0.338 
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Allow  for  a  Finite  Diode: 


Displacement  Vector: 


z-component  of  the  Electric  Field: 


25 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


26 


LIST  OF  REFERENCES 


1.  Anderson,  Weston.  "Role  of  Space  Charge  in  field 

emission  cathodes".  Journal  of  Vacuum  Science 

Technologies.  Volume  12  Number  2.  Mar/Apr  1993. 

2.  Barker,  Robert  and  Schamiloglu,  Edl .  High-Power 

Microwave  Sources  and  Technologies.  IEEE  Press,  New 
York.  2001 


3.  Griffiths,  David.  Introduction  to  Electrodynamics. 
Third  Edition.  Prentice  Hall,  Inc.  New  Jersey.  1999. 

4.  Goldsten  and  Rutherford.  Introduction  to  Plasma 
Physics .  IOP  Publishing,  1995. 


5. 


Jackson,  John  David.  Classical  Electrodynamics. 

Third  Ed.  John  Wiley  and  Sons,  New  Jersey,  1999. 


6.  Kirstein,  Kino,  and  Waters. 
McGraw-Hill,  New  York,  1967. 


Space  Charge  Flow. 


7.  Langmuir,  Irving  and  Blodgett,  Katherine.  "Currents 
Limited  by  Space  Charge  between  Concentric  Spheres". 
Physical  Review  Letters .  Volume  24,  1924. 

8.  Lau,  Y.Y.  "Simple  Theory  for  the  Two-Dimensional 
Child-Langmuir  Law".  Physical  Review  Letters .  Volume 
87,  Number  27.  December  2001. 


9.  Luginsland,  John  et  al .  "Beyond  the  Child-Langmuir 
law:  A  review  of  recent  results  on  multidimensional 
space-charge-limited  flow".  Physics  of  Plasmas,  Vol . 
9,  Number  9,  May  2002. 

10.  Rokhlenko,  A.  and  Lebowitz,  J.  "Space-Charge-Limited 
2D  Electron  Flow  between  Two  Flat  Electrodes  in  a 
Strong  Magnetic  Field".  Physical  Review  Letters . 
Volume  91,  Number  8.  August  2003. 


27 


11.  Umstattd,  Ryan  and  Luginsland,  John. 
Space-Charge-Limited  Emission: 

Characteristics  and  Applications". 
Letters.  Volume  87,  Number  14.  October 


"Two-Dimensional 
Beam-Edge 
Physical  Review 
2001 . 


28 


INITIAL  DISTRIBUTION  LIST 


1.  Defense  Technical  Information  Center 
Ft.  Belvoir,  Virginia 

2 .  Dudley  Knox  Library 
Naval  Postgraduate  School 
Monterey,  California 

3.  Captain  Ryan  Umstattd,  USAF. 

Department  of  Physics 
Monterey,  California 

4.  Professor  Chris  Frenzen 
Department  of  Mathematics 
Monterey,  California 

5.  Chairman,  Physics  Department,  Code  PHMW 
Naval  Postgraduate  School 

Monterey,  California 

6.  Engineering  and  Technology  Curriculum  Office,  Code  34 
Naval  Postgraduate  School, 

Monterey,  California 


29 


